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We study the thermodynamics of the diffusive SFS trilayer composed of thin superconductor (S) 
and ferromagnet (F) layers. On the basis the self-consistent solutions of nonlinear Usadel equations 
in the F and S layers we obtain the Ginzburg-Landau expansion and compute the condensation free 
energy and entropy of the 0 (even) and tt (odd) order parameter configurations. The first order 0 —tt 
transition as a function of temperature T occurs, which is responsible for a jump of the averaged 
magnetic field penetration depth \{T) recently observed on experiments [N.Pompeo, et. al., Phys. 

Rev. B 90, 064510 (2014)]. The generalized Ginzburg-Landau functional was proposed to describe 
SFS trilayer for arbitrary phase difference between the superconducting order parameters in the S 
layers. The temperature dependence of the SFS Josephson junction critical current demonstrates 
the strong anharmonicity of the corresponding current-phase relation in the vicinity of the 0 — tt 
transition. In rf SQUID, coexistence of stable and metastable 0 and tt states provides integer and 
half-integer fluxoid configurations. 

PACS numbers: 74.45.-l-c, 74.78.Na, 74.78.-w 


I. INTRODUCTION 

The ground state of the superconductor-ferromagnet- 
superconductor (SFS) trilayer at zero current can be 0 
or TT state, depending on the value of the phase differ¬ 
ence between the superconducting order parameters in 
the two S electrodes. This phenomenon is related to 
the damped oscillatory behavior of the Cooper pair wave 
function in the ferromagnet due to the proximity effect^i^ 
(for more references and reviews, see Refs. ii). Usually 
experiments directed towards the observation of the 0 — tt 
crossover in SFS trilayer were concentrated on the mea¬ 
surements of the critical Josephson current Ic of the SFS 
junction^5^. The 0 — tt transition manifests itself in the 
vanishing of Ic if higher-order harmonics of the current- 
phase relation are negligibleiSrJ^. 

Recently, an unusual electromagnetic response of SF 
systems was reported as a manifestation of the Cooper 
pair wavefunction oscillations inside the ferromagnet. 
Measurements of the London penetration depth in thin 
Nb/Ni bilayersi^ reveal a slightly nonmonotonic depen¬ 
dence of the penetration depth on the F layer thickness, 
which was in accordance with the theoretical analysis^^. 
Anomalous Meissner effect in hybrid SF structures was 
the subject of several theoretical worksi^^— predicting an 
unusual paramagnetic response of such systems. Vanish¬ 
ing or inversion of the Meissner effect is believed to be at¬ 
tributed to spin-triplet superconducting correlation o^^d^ 
generated in inhomogeneous F layer due to proximity 
effect and should result in the in-plane Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) instability22i“— . Unusual 
drop of the screening with decrease of temperature was 
observed recently in22 by microwave measurements of the 


London penetration depth A in Nb/Pdo. 84 Nio.i 6 /Nb tri¬ 
layers. 

The transition temperature of SF structures into the 
normal state has been examined both theoretically^lvSi 
and experimentally^^— (see Ref. 1^1^. for reviews). 
However, the study of the thermodynamic properties of 
the phase transition between 0 and tt states of SF hybrids 
is more sparse. A first-order 0 — tt transition was pre¬ 
dicted for diffusive SFS junctions with a homogeneous F 
barrier, if the current-phase relation takes into account 
the second harmonic contributiooii. Experimental ev¬ 
idence of a 0 — TT transition in SFS (Nb/CuxNii_x/Nb) 
Josephson junction was obtained from the measurements 
of the tenmerature dependence of the critical current 
in Refs. Note that a small modulation of the 

thickness of the barrier may favor the continuous 0 — tt 
transitionii. A first-order transitions between 0 and tt 
states by temperature variation were demonstrated in 
both the clean^iiS^ and dirty^l limits using numerical 
self-consistent solutions of the microscopic Bogoliubov- 
de-Gennes^i22 or Usadel^ equations, respectively. The 
0 — TT transition in ballistic SF systems has been shown 
to have a pronounced effect on the distribution of the 
Cooper pair wavefunction in the F region and the am¬ 
plitude of order parameter A(T) in S-layer s^"^i^^ . The 
theoretical model Proposed in Ref. argued that the 
observed jump A(T) in Nb/Pdo. 84 Nio.i 6 /Nb trilayers is 
related to the first order phase transition from 0 to tt 
state on cooling. 

In this work we develop a theoretical approach based 
on the nonlinear Usadel equations providing a general 
description of diffusive SFS junction with thin supercon¬ 
ducting layers at the transition from 0 to tt state. The 
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leakage of the Cooper pairs weakens the superconductiv¬ 
ity near the interface with a F metal due to the proximity 
effect^. The magnitude of the superconducting order pa¬ 
rameter suppression depends on the parameters charac¬ 
terizing the system such as the SF interface transparency, 
the thickness of the S and F layers, etc. For large inter¬ 
face transparency this effect seems to be especially strong 
and results in suppression of the superconducting order 
parameter and the transition temperature Tc of a thin su¬ 
perconducting layer in contact with a ferromagnet metal. 
If the thickness of a superconducting layer is smaller than 
a critical one, the proximity effect totally destroys a su¬ 
perconductivity. From the self-consistent solutions we 
obtain the Ginzburg-Landau expansion and compute the 
condensation free energy and entropy of the possible or¬ 
der parameter configurations as a function of tempera¬ 
ture T. As T varies, we find that the first order phase 
transition between 0 and tt states occurs, which is respon¬ 
sible for a jump of the averaged penetration depth X{T) 
observed in Ref. [2^. We also calculate the current-phase 
relation I (ip) of the SFS junction which reveals strong 
contribution of the higher harmonic terms. The 0 — tt 
states coexistence and switching leads to new modes of 
magnetic flux penetration in superconducting loop con¬ 
taining the SFS junction. 

The paper is organized as follows. In Sec.lUwe briefly 
discuss the basic equations. We analyze the case of thin 
S layers and obtain the approximate solutions of the non¬ 
linear Usadel equations in F layer near the superconduct¬ 
ing critical temperature T < Tc- In Sec. IIIII we find the 
temperature of the second-order superconducting 
phase transition to the normal state for 0 and tt order- 
parameter configurations. The Sec. US is devoted to the 
analysis of the temperature-driven transition between 0 
and TT states. We obtain the Ginzburg-Landau expan¬ 
sion and find the ground states of SFS junction near the 
critical temperatures In Sec. Ewe generalize the 

Ginzburg-Landau description for arbitrary phase differ¬ 
ence If between the superconducting order parameters of 
the S layers and find the strongly nonsinusoidal current- 
phase relation of the SFS junction in the vicinity of 0 — tt 
transition. In Sec. ED we show that the coexistence of 0 
and TT states leads to peculiarities of the magnetic flux 
penetration in superconducting loop with a single SFS 
junction. Sec. IVII| contains a brief summary and discus¬ 
sion. 


II. MODEL AND BASIC EQUATIONS 

Let us consider a SFS trilayers with a transparent SF 
interfaces and thin S layers of thickness dg ^ where 
is the superconducting coherence length. The con¬ 
sidered geometry of the SFS structure is presented in 
Fig. [TJ In the previous observations of O-tt transition 
in SFS junctions by temperature variation^ the super¬ 
conducting electrodes were rather thick to overcome the 
pair-breaking proximity effect of F layer and then the 



FIG. 1: (Color online) The schematic behavior of the pair 
wave function F{x) = Fsj inside the SFS trilayers. The red 
solid line represents approximately the behavior of the pair 
wave function in an even mode (0-state). Due to symmetry 
the derivative dxFf is zero at the center of F layer. The pair 
wave function in the odd mode (blue dashed line) vanishes at 
the center of F layer, and F(x) has a tt— shift in diametrically 
opposite points (tt— state). 


influence of the O-tt transition on the superconducting 
order parameter A(T) in the electrodes was negligible. 
From the theoretical point of view, weak depairing in S 
electrodes means that the pair potential at the SF in¬ 
terfaces is equal to its bulk value (the so-called rigid¬ 
boundary condition^). Here we study the SFS structure 
with relatively thin S layers and demonstrate that the O-tt 
transition leads to the jump of the amplitude of the su¬ 
perconducting order parameter, providing the anomalous 
temperature behavior of the effective penetration depth 
of the whole structure. 

To elucidate our results we start the qualitative discus¬ 
sion of the proximity effect on the properties of SFS sand¬ 
wiches if the thickness of S layers is small enough. Due 
to the damped oscillations of the pair wave function F 
in the ferromagnetic layer, two different order-parameter 
configurations are possible inside the SFS trilayers (see 
Fig. HI). The first one corresponds to the case than the 
order parameter is an even function of the coordinate x, 
chosen perpendicular to the layers, and does not change 
its sign in F layer. It means that in the ground state the 
superconducting phase in both S layers must be the same 
(0—phase). For the second case the pair wave function 
is odd in x and cross zero at the center of the F layer 
which causes a tt— shift in the superconducting phase of 
different S layers. This configuration corresponds to the 
TT— phase of the SFS structure. Since pair-breaking prox¬ 
imity effect depends on the structure of the Gopper pairs 
wavefunction in ferromagnetic layer, a suppression of su¬ 
perconductivity in thin S layers is expected to be dif¬ 
ferent for 0 and tt states. As a result the equilibrium 
superconducting gaps Aq^tt for 0 and tt states are dif¬ 
ferent (Ao(T) ^ A^(T)), and this leads to the different 
effective penetration depth A(T) ~ 1/A(T) in 0 and tt 
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states. So, O-tt transition in SFS trilayers has to be ac¬ 
companied by a jump of A. The coexistence of stable and 
metastable states in the vicinity of 0 — tt transition leads 
to a strong anharmonicity of the current-phase relation 
in the SFS junction under consideration. As a result, 
peculiarities of the magnetic flux penetration in a meso¬ 
scopic superconducting loop containing the SFS junction 
are expected. 

Taking in mind that the F interlayer is dilute ferromag¬ 
netic alloys like CuxNii_x or PdxNii_x we use the Usadel 
equations^^ which are convenient in the diffusive limit 
(see^ for details). Moreover it is important to take into 
account the magnetic disorder which is already present 
in the magnetic alloys and provides the main mechanism 
of the temperature induced O-tt transitionSi^. We de¬ 
scribe it by introducing the magnetic scattering time Tg. 
Note that in practice the exchange field h acting on the 
electron’s spins in the ferromagnet and magnetic scatter¬ 
ing rate are much larger than the superconducting 
critical temperature T^- h, ^ Tc. The complete 
nonlinear Usadel equations for the normal G{x, w, h) and 
anomalous F{x, uj, h) Green’s functions in F layer ar o^'^^ : 


Df 


[G{x, uj, h) d‘^F{x, UJ, h) — F{x, uj, h) d^G(x, uj, h)] 


+ 


ih - 


G{x, UJ, h) 


F{x, uj,h) = 0 , (1) 


© for the S layers have to be completed by the self- 
consistency equation for the superconducting order pa¬ 
rameter A(a;): 


A{x) ='!rTp'^Fs{x,uj), (6) 

UJ 

where p is BCS coupling constant and T^q is the critical 
temperature of a bulk sample of the material S. 

For 0 and tt states of SFS trilayers we have 
F'^{x,uj,h) = F*{x,uj,—h) = F(x,uj,h), and one can 
replace F^{x, uj, h) by F{x, uj, h) in Eqn. ([2]). Just below 
the critical temperature anomalous Green’s functions 
are small and the condition ([5]) can be rewritten as 

G{x,uj,h) — F‘^{x,uj,h)/2 . (7) 

For ds is, where is = \JFJsI^.-kTcui is the supercon¬ 
ducting coherence length, the variations of the functions 
9s{x) and A(x) in the superconducting layers are small: 
9s{x) ~ 9s, A(x) ~ A. Therefore, we can average Eq. ([3]) 
over the thickness of the S layers, using the boundary 
condition ([S]). Finally, we obtain the following boundary 
condition: 


d9s 

dsif 

/ A cos 9s — UJ sin 9s A 

ds 

SJ Ss 

V ^Tco j 


G^{x,uj, h) + F{x,uj, h)F^{x,uj, h) = 1. (2) 


ds 


-s/ 





(9) 


Here Df is the diffusion constant in the ferromagnet, 
UJ = 27rT(n -I- 1/2) is a Matsubara frequency at the tem¬ 
perature T. The equation for the function F'^(x,uj, h) = 
F*{x,uj,—h) coincides with Eqn. ([T]).^. For the 0—state 
we should choose the even anomalous Green’s functions 
F while for the tt— state it should be the odd one. Using 
the usual parametrization of the normal and anomalous 
Green functions Gg = cos 9s and Fg = sin 9s, the com¬ 
plete nonlinear Usadel equation in superconucting layers 
can be written for w > 0 as^ 

— ^^d‘l9s-\-ujsin9s = Acos9b , (3) 

where Dg is the diffusion constant in the superconductor. 
Assuming the SF interfaces to be transparent we have at 
X = ±df/‘M: 


Fg = F, cjg d^Fg = af d^F , (4) 

where cr/ and Us are the normal-state conductivities of 
the F and S metals, respectively. The boundary condition 
at the outer surfaces x = ±(c?s -I- df/2) is 

ds,Fs = 0 . (5) 


For thin S-layers ds is = \/dJJ2ttT^, the inverse 
proximity effect is substantial, and the Usadel equation 


where s = x/if and s/ = df/2if. The top (bottom) 
sign in corresponds to the 0—phase (tt— phase), re¬ 
spectively. 

Applying the methodic we can find the approximate 
solution of the non-linear Usadel equation © in F layer 
near the superconducting critical temperature T < T^. 
For the 0—state we should choose the even anomalous 
Green’s functions F while for the tt— state it should be 
the odd one: 


F{s,uj) 


a cosh(gs) 
bsi'Dh{qs) 


( 10 ) 


1 / 3A j cosh(3fcs), 0 — phase 

8k‘^ \ 4 J I 5^ sinh(3fcs), tt — phase 


where 0 = Df/h, = 2(yjlh + i + a sgn(w)) and a = 
I/ts/i is the dimensionless magnetic scattering rate. The 
complex wave vector q is determined by the relations: 

=F (a + V4) (H) 

for 0 and tt phases, respectively. If Tc < h, we 
may neglect the Matsubara frequencies in the Eq. o 
assuming that = {ki + ik 2 )‘^ = 2(q; -I- f) for w > 0 : 


ki 




\/l + a'^ + a. 



( 12 ) 
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Then the decay characteristic length and the oscilla¬ 
tion period ^f 2 , rnay be written as 

5/2= 0/^2- (13) 

The ratio of the characteristic lengths ^/i /^/2 < 1 
clearly shows that magnetic scattering decreases the de¬ 
cay length and increases the oscillation period^. 


A. Even mode (0—phase) 


In the limit |a| <C 1 we obtain from the even solution 
(I10lll|) the expansion of F in powers of the amplitude a: 

F{s) ~ fcosh{ks)-fgo{s)+o{f^), (14) 

go{s) = -^[4:ks{a + i/A)smh{ks) 

-|-(a-I-3i/4) cosh(3fcs)] . 

Using the first boundary condition ([4]) at s = s/ we get 
the relation between the amplitude a and Green’s func¬ 
tion Fs = sin 0s in superconductor: 


/ 


fo + f^ 


9o{sf) 
cosh(fcs/) ’ 


/o 


F, 

cosh(fcs/) 


(15) 


Substitution of the solution o to the relations dm) 
results in the following equation with respect to the am¬ 
plitude Fg in the S layers: 

(cc + 1/ro) = A - i {^F! + sAoF!) (16) 

where 


Tq ^ = eirTcok tanh(/cs/) (17) 

is the depairing parameter of even mode and 


Aq 


t^Tco 


^tanh(fcs/) 


a + z/4 tanh(fcs/) 
k cosh^(fcs/) 


Sfja + i/A) 
cosla^ (ksf) 


(18) 


Here the key parameter 


e = 


CTs dg^f 


determines the influence of the proximity effect on the S 
layers. 



4 A, 


FIG. 2: (Color online) The typical dependence of the criti¬ 
cal temperature T^’’" on the thickness of F layer dj for even 
mode(0—phase) (solid red line) and for odd mode(7r—phase) 
(dashed blue line). Here we choose: dg = 2^s; = 0.12; 

is/is = 3 (e = 0.18); hrg = 7. 


Using the first boundary condition (|T]) at s = s/ we get 
the relation between the amplitude b and Green’s func¬ 
tion Fg = sin 9g in superconductor: 


f = u- f, 


3 

sinh(A:sj) 


, U = 


Fg 


sinh(fcs/) 


( 20 ) 


Substitution of the solution m to the relations (I5P) 
results in: 

(uj + 1/rO = A - i (AF," + eA^F!) , (21) 


where 


^ = enTcok coth(fcs/) 

is the depairing parameter of odd mode and 


( 22 ) 


— coth(A:s/)-h 


a-|-j/4 coth(fcs/) Sf{a + i/4:) 


k sinh^(fcs/) sinh"‘(fcs/) _ 


(23) 


B. Odd mode (tt— phase) 

In the limit |5| <C 1 we obtain from the odd solution 
(fTOTTTt the following expansion of F in powers of the 
amplitude b: 

F{s) ~ fsinh{ks)+fg^{s) + o{f^), (19) 

S 7 r(s) = ^ [4A:s(a-b f/4) cosh(fcs) 

— (a-b 3f/4) sinh(3fcs)] . 


III. THE CRITICAL TEMPERATURE OF SFS 
TRILAYER 


To find the temperature of the second-order su¬ 
perconducting phase transition, the equations dm), (EiD 
should be linearized with respect to the Fg 1 


Fg ~ Fgo — 


A 

W -b 


(24) 
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Substituting Eq. (|24l) into the self-consistency equation 
m one obtains the equation for the critical temperature 

Y^0,7r. 

c 


In 




T, 


cO 




(25) 


where 'I' is the digamma function. The depairing param¬ 
eter = l/27rTTo,^ 


J ^ tanh(A:s/), 0 - phase , , 

2 (A: coth(A:s/), tt — phase 

is responsible for the superconductivity destruction in the 
S layers due to the proximity effect. Figure [2] shows a typ¬ 
ical dependency of the critical temperature on the 
thickness of F layer df, obtained from Eqs. (1^ . (1^ . 
The crossing of the curves T^{df) and T^{df) occurs 
at dp and for df > d*j the critical temperature of the 
TT— phase becomes higher than the critical temperature 
of the 0—phase. 





FIG. 3: (Color online) The typical dependence of the coef¬ 
ficients of the Ginzburg-Landau expansion aP’’^, h^’’^ on the 
thickness of F layer df. (6°) - closed (open) red circles; 

{h") - closed (open) blue triangles. Here we choose the 
parameters for Fig. (2] 


IV. PHASE TRANSITIONS IN SFS TRILAYERS 


To obtain the Ginzburg-Landau (GL) expansion for 
the 0— and tt— states near the critical temperature 
we can use the equations (ITO1) . (I^ . Assuming \Fs\ <C 
1 and using Fgo (IMl) as a zero-order approximation we 
find the solution of Eqs. (HU),® within the first-order 
perturbation theory: 


Fs = sin Os ^ Fso 


2 _ LU + \/ 


(27) 


Substitution of (l27)) into the self-consistency equation ([6]) 
one obtains a dependence of the superconducting gap A 
on the temperature T = T^’’" — ST 


+ = 0 , 


T 


(28) 


where the coefficients and 5°’'^ are determined by the 
following expressions: 


0,7r _ 


= 


1 — Re 
-1 


L!o,,TW(l/2 + r!o..) 


i0,ir\2 


Re 


T(2)(l/2 + fIo..) 


(47rr; 

DTT 1 c 


(29) 

(30) 


where 'I'*(")( 2 ;) = d^'^>{z)/dz’^ . Natnrally all parameters 
jiO.TT (different for 0 or tt states. Fig¬ 
ure [3] shows a typical dependence of the coefficients of 
the Ginzburg-Landau expansion on the thick¬ 

ness of F layer df, obtained from Eqs. ®,(|3ni)- 


The equilibrium value of superconducting gap 
, n 0°’’^ ry -T 

^O.TT - yO.TT (31) 

corresponds to the extremum of the standard Ginzburg- 
Landau functional 


F^'l{T)=Eo 


rp _ 


^O.TT 


T 


0,7r 




(32) 


where the characteristic energy Eq = N{0)SjdT^Q is de¬ 
termined by the total electron density of states N{0), the 
cross section area Sj of the junction, the total thickness 
of the trilayer d = 2ds + df and the critical temperature 
Tcq. The value of the superconducting order parameter 
A, the temperatures T, T^’’^ are assumed to be measured 
in the units of Tcq. The functional (15^ provides us the 
complete description of the ground states of SFS junc¬ 
tion near the critical temperature . The equilibrium 
energy of the system E°’'^{T) = FQ^(Ao, 7 r) is 


E°'^{T) = -Ac 


[a°’^ {T°’^ - T)/T°- 
2b°’^ 


(33) 


The condition of the first order transition between 0 and 
TT states is F°(Ao) = F'^(A^), and then the temperature 
of this transition Tq is determined by: 


T° -To _a^ [¥ 
T^ - To ~ ¥ 


(34) 


df ^ d*p first the transition from normal state to 
the superconducting 0 state occurs, but further decrease 
of the temperature provokes the transition from 0 to tt 
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FIG. 4: The (T, df) phase diagram for the SFS trilayers in the vicinity of crossing (d/ « df) of the curves T^{df) and (df): 
(a) d*f ~ 1.95^/; (b) d*f ~ 5.31^/ . At temperature To the first-order transition between the 0— and tt— states takes place. 
Here we choose the parameters for Fig. [SjThe insert gives the dependence of the value of the superconducting gap jump 
Ao(ro)/A^(ro) (1371) and the latent heat Q (1361) on the transition temperature To- 


state. The correspondent (T, df) phase diagram for the 
SFS trilayers is shown in Fig. 0] 

The 0 — TT transition is accompanied by a discontinu¬ 
ity in the entropy S^’'^{T) = — [9if°’'^(T)/9T] at the 
temperature Tq: 

S^Tq) ^ - To 

S°{To) T^-To' ^ ’ 

Then the latent heat at the first order 0 — tt transition is 


Q = ±To [5^(To) - ^°(To)] > 0 (36) 

for 0 TT transitions, respectively. Simultaneously, at 
the transition temperature Tq the superconducting order 
parameter jumps from Aq to Att or vice versa. The ratio 
of values Aq and Ajr is given by the expression: 


A^(To) 


a- /T--To\ 

6- V ) 



(37) 


The inserts in Fig. 0] show the dependence of the ra¬ 
tio Ag(To)/A^(To) (1571) and the latent heat Q{Tq) (1551) 
on the transition temperature Tp. Certainly, the jump 
of the superconducting gap provokes the jump of the 
London penetration depth 1/A ^ A. Figure [S] shows 
schematically the temperature dependence of the equi¬ 
librium gap Aq ^ (1311) and the Ginzburg-Landau energy 
(1551) in the vicinity of the 0 —>■ tt transition for 
the case T/’ < T° and > W. We readily see that the 
superconducting order parameter decrease results in the 
positive jump of the London penetration depth. At all 
reasonable parameters h, we obtain namely this 

scenario, which is indeed realized on the experiment^. 


V. THE CURRENT-PHASE RELATION 



FIG. 5: (Golor online) Schematic dependence of the gap Aq 
(EB and the Ginzburg-Landau energy FGL{Ao^,r) (l33l) on the 
temperature T: < T® and [a?/h)'” > /b)^ . The inset 

gives the schematic temperature dependence of the penetra¬ 
tion depth A ~ I/Aq.tt. 


from 0 to TT state. The theory describing how the 0 state 
is transformed into the tt state for diffusive SFS junctions 
was developed in Ref. [ij for rigid boundary conditions^ 
at SF interface. It was shown that the critical current 
does not vanish at the transition and is determined by 
the second-harmonic term in the current-phase relation 


Let us now discuss the peculiarities of the Josephson 
effect in the SFS trilayers at the first-order transition 


I (ip) = Ii sin p + I 2 sin 2p , 


(38) 
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which corresponds to the following phase-dependent con¬ 
tribution to the free energy of the junction: 

h 

—Ii (1 — cos (fi) —;7 (1 ~ cos 2ip) , 

(39) 

where $o = 7r/ic/e is the flux quantum. Since the am¬ 
plitudes of the harmonics Ii and I 2 depend strongly on 
the superconducting order parameter A in S layers (see 
Appendix for details), the current-phase relation of 
the junction I (if) seems to be very sensitive both to a 
suppression of superconductivity in thin S layers and to 
jumps of A during the transitions between 0 and tt states 
of the SFS trilayers. 

To consider the Josephson coupling in the SFS trilay¬ 
ers in the presence of first-order 0 — tt transition, the 
Ginzburg-Landau functional (|32|) should be generalized 
for arbitrary phase difference (p between the supercon¬ 
ducting order parameters in the S layers. For two coupled 
S layers with the order parameter Ai ^2 = Aexp{±iip/2) 
the Ginzburg-Landau expansion Fqj^ of the free energy 
includes the mixing quadratic term A 1 A 2 + A*A 2 ^ 
A^ cos if. The mixing terms of fourth order have the 
form (|Aip + |A 2 p)(AiA 2 + A);A 2 ) ^ A'^cos:/? and 
A^(A 2 )^ -I- (A*)^A 2 ^ A^cos2(/j, and provide both ip 
and 2i^ periodicity of the function Fqj^. Therefore, in 
general, the Ginzburg-Landau expansion can be written 
as 


Ej{p) = 

2ttc 


F^^{T, p)/Eo = -a^{T)A^ + ^A^ , (40) 

which incorporates the phase-dependent contribution via 
the coefficients 

Tc — T 

=71 cos(/J-f "" (1-b72COsyi), 

^ C 

= Po + Pi cos p + P2 cos(2(/j). 

For convenience the superconducting order parameter A, 
temperatures T, Tc, T* and the exchange field h are mea¬ 
sured in units of Tcq. To find the parameters Tc, T*, 71^2, 
/3 o,i, 2 we take into account that the generalized functional 
(HOI) reduces to the expression (1001) if (/; = 0, tt: 

F^^{T, 0) = FO^(T) , tt) = F5^(T). (41) 

As a result for T < min(T° ,T^) the coefficients Tc, T*, 
71 . 2 , can be expressed as 


/yi* _ 


a°/T° + a-^/T^ 


rp'* 

Tc=^ (a° 
2 ^ 




+ a'^) , 


71 = 




72 D \ 2^o 2^-77 


Pi = ib° - b^) 12, p^+p2 = ( 6 ° + 5 -) /2 . 


The coefficient P 2 in expansion (HOI) can be determined 
from the current-phase relation for rigid boundary con¬ 


ditions (see Appendix lAT): 


P2 = -TTP; 


h C/. 

192 T3 d 


1 


S .a 
2 sinh 6 


*/4 


k sinh^ 6 
cosh 6 — 


sinh^ 


(42) 


where 6 = 2sfk = d/Z^/i + idf /^f 2 - Substitution of the 
equilibrium value of the order parameter 

Al{T) = a^{T)/b^. (43) 

into the expression (j40p provides the temperature depen¬ 
dence of the free energy E{T, p) of the SFS trilayers for 
an arbitrary phase difference p 

E{T,p)lE^ = -al{T)l2b^. (44) 

which results in the following current-phase relation 

I{p) = {2elh)dEldp: 

I{T,p) = Jo sin p X (45) 

71 + 72 ^ y^ ^ {Pi + 4/32 cos p) , 

where Jq = 2ttcEq/^q. The current-phase relation I{p) 
(H51) is shown in Fig. for fixed thickness of the bar¬ 
rier df = 1.94^^- in the vicinity of the 0 — tt transition 
and several values of the temperature T. For chosen pa¬ 
rameters of SFS trilayers (see Fig. [3|) the first O-tt tran¬ 
sition takes place at df < 1.95^/. Corresponding value 
of the coefficient P 2 — —0.013 gives a rather large value 
of the second-harmonic term (I/ 2 I ^ |Ji|)j which domi¬ 
nates near the transition. Due to the strong contribution 
of higher harmonics the current-phase relation is rather 
anharmonic, and the maximal value of |J(T, p)\ occurs 
at p TT 12. Figure [Bh shows the temperature depen¬ 
dence of the critical current Jc(T) = maX(^ \I{T, p)\ near 
the temperature Tq- The characteristic multimode an- 
harmonicity of the current-phase relation results in the 
disappearance of the typical nonmonotonic temperature 
dependence of the critical current in a vicinity of the 0— tt 
transition. We have obtained the positive amplitude of 
second harmonic I 2 > 0, which means that it occurs dis- 
continuously by a jump between 0 and tt phase states at 
the transition point Tj, where the critical current Ic(T) 
formally changes it sign. The shift of the temperature 
Tj with respect to Tq depends on the higher harmonics 
contribution. 

For comparison figure 0 shows the current-phase rela¬ 
tions I{p) (l45l) and the temperature dependence of the 
critical current Ic{T) in a vicinity of the second tt— 0 tran¬ 
sition (see Fig. El which takes place at df < 5.31^/. Gor- 
responding value of the coefficient P 2 — —0.18 10“"^ gives 
a small value of the second-harmonic term (I/ 2 I <C |Ji|), 
and the the current-phase relation is quite harmonic ex¬ 
cept the case T Ri Tq ~ T/. The SFS junction reveals the 
typical nonmonotonic behavior of the critical current Ic 
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FIG. 6: (Color online) (a) Current-phase relation of SFS junction I{ip) (I45II for several values of the temperature T/Tco = 
0.30, 0.31, 0.32 in the vicinity of the 0 — tt transition for df = 1.94^/ {T°/Tco — 0.332, jTco — 0.327, T^jTca — 0.329). The 
case T = Tq ~ 0.311 Tco is shown by the dash-dotted curve (/a = 7o A§(T), where Ao(r) = 1.764 Tc tanh(1.74^rc/r — 1 is 
BCS superconducting gap for the temperature T). (b) Dependence of the critical current Ic{T) = max|7(T, (p)| (red solid line) 
and the amplitudes of harmonics 7i, I 2 (blue dashed lines) on temperature T. The insert gives the temperature dependence 
of a relative amplitude of the first harmonic k = Ii/l 2 - Here we choose the parameters of Fig. [3] ds = 2^s; ojjoa = 0.12; 

= 3 (d} ~ 1.946C/, £ = 0.18) and h = lOTco, Kts = 7. 




(p / n 


FIG. 7: (Color online) (a) Current-phase relation of SFS junction 7(i/j) (I45II for several values of the temperature T/T^o = 
0.28, 0.29, 0.30 in the vicinity of the second 0 —tt transition for d/ = 5.3^/ (T ®/Tco — 0.317, /Tco — 0.318, Tc/Tco — 0.3175 ). 
The case T = To — 0.29 Tco is shown by the dash-dotted curve (7a = lo Ao(r), where Ao(T) = 1.764 Tc ta,Ti\i{1.7A^jTc/T — 1 
is BCS superconducting gap for the temperature T). (b) Dependence of the critical current Ic{T) = max|7(r, (p)| (red solid 
line) and the amplitudes of harmonics 7i, I 2 (blue dashed lines) on temperature T. Here we choose the parameters of Fig. |3] 
da = 2^,; af/aa = 0.12; ^a/0 = 3 (d} ~ 5.31^/, £ = 0.18) and h = 10T,o, liTa = 7. 


as a function of the temperature T, and the position of 
the cusp Tf naturally coincides with the temperature Tq. 

If we restrict our consideration to the two harmonics 
approach (13811 . one will be able to determine the ampli¬ 
tude of both harmonics 7i and I 2 via critical current 7c 
measurements, as it has been proposed in Ref. |4l|. Since 
a relative amplitude of the first harmonic k = I 1 /I 2 is 
small (|k| < 2) in the vicinity of 0 — tt transition (see the 
insert in Fig.|B]D), the system has two stable states ip = Q 
and (/? = TT at 7 = 0. To ’’depin” the Josephson phase 


from the low energy 0 (tt) state or from the high energy 
TT (0) for 7i > 0 (7i < 0), respectively, the critical current 

/c±(aj) = ^ (v^K2 + 32±3|K|j (46) 

X ^\/-(- 32 =p |k|^ 

should be applied^. During switching from the voltage 
state to the zero resistance state at T Ri Tq, the phase 
may stick in the high energy state with a probability close 
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FIG. 8: (Color online) The shape of the Josephson vortex 
V5o,7r(a;) (14811 . (1491 1 for two values of the relative amplitude of 
the first harmonic k = I 1 /I 2 in the vicinity of the 0 — tt 
transition: blue dashed line - |«:| = 0.05; red solid line - 
|k| = 10“^. 

to 50%, and the smaller current Ic- becomes measurable 
in the SFS structure with a low damping^. Then for 
|«:| <C 1 near 0 — tt transition the amplitude of the har¬ 
monics /i ,2 is determined by the relations 


FIG. 9: (Color online) Magnetic flux through the single¬ 
junction loop $ as a function of the external flux 4>e for the 
normalized inductance Li = 250 and several values of the 
temperature T/Tco = 0.28, 0.29, 0.3, 0.31, 0.32 in the vicin¬ 
ity of the 0 — TT transition for df = 1.94^/. The case for 
T = To — 0.31 Tco is shown by the dash-dotted curve. Here 
we choose the parameters of Fig. [6] and Li — 500. 

ground state of the circuit is determined by minimizing 
the SQUID free energy 


|/iI ~ (4+ -Ic-)/V2, h- (/c+ + Ic-) /2 . (47) 


The structure of the Josephson vortex which may ex¬ 
ist in a long Josephson junction with a large second har¬ 
monic (|k| < 2) is rather peculiar and is described by the 
following expressionsii 


To = 


2tt - TKix/Xj2), X <0 
Tk{x/Xj2), x>0 


, Ttt =^0-^^ (48) 


for /i > 0 (/i < 0), respectively, where 

2(1 -I- |«;|) 


Tk{x) = arccos 1 - 


1 -h |K|cosh^(x\/l -h |k|) 


(49) 


and the Josephson length Xj 2 = c^oS/Sn'^ht depends 
on the current density I 2 /S and the effective junction 
thickness t. The change of the form of the Josephson 
vortex in the vicinity of the 0 — tt transition is shown in 
Fig. m 


VI. SINGLE-JUNCTION LOOP 


W{T, cl,) = E{T, </>) + ^ (<(, - ^,)2 , (50) 

where (^e = 27r$e/<i)o - is the normalized magnetic flux of 
an external field through the loop^. The total magnetic 
flux through the loop $ = (0/27r) $0 is related to the 
external flux as 


^e=^ + LiI{T, ct))/h , (51) 

with the normalized inductance Lj = 2ttL Iq/ 
Figure [S] shows the dependence 4)(<i)e) determined by 
Eqs. (HOD and dSU for several values of the tempera¬ 
ture T near the transition temperature Tq. New features 
of 4)(<i)e) dependence appear for T close to Tq due to 
the coexistence of stable and metastable 0 and tt states. 
Strong anharmonicity of the current-phase relation for 
T Ki Tq Tc results in the coexistence of integer and 
half-integer fluxoid configuration in SQUID’s and gener¬ 
ation of two flux jumps per one external flux quantum. 
This behavior is similar to the magnetic-flux penetration 
in superconducting loop with a clean SFS junction at low 
temperatures T <C T^. 


Let us consider a small superconducting loop with an 
inductance L interrupted by the single SFS junction (rf 
SQUID). We assume that the junction is described by 
the current-phase relation (H51) and choose the parame¬ 
ters in a narrow region near the 0 — tt transition. The 


VII. SUMMARY 

We have studied the thermodynamics of diffusive 
SFS trilayer with relatively thin S layers through self- 
consistent solutions of nonlinear Usadel equations, in 
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the dirty limit. Our results may be viewed as gen¬ 
eralization of those obtained in Ref. when the su¬ 
perconducting electrodes are the rather thin, and the 
critical temperature is affected by the ferromagnetic 
layer. We have shown that as the temperature T is 
varied a given SFS junction can flip from the 0 state 
to the TT state. The resulting phase transition is first- 
order, in agreement with the experiments^i^i^ and is 
responsible for a jump of the amplitude of the super¬ 
conducting order parameter A, providing the anomalous 
temperature behavior of the effective penetration depth 
A(T) observed in Ref. [H. Taking the typical parameters 
N{0) ^ 10^^ erg“^ sm“^, Tco ^ lOK Ri 10“^^ erg we may 
estimate for the Nb/Pdo. 84 Nio .le/Nb trilayer.2M5 ~ 
30 nm, S' ^ 0.1 sm^) the latent heat Q ^ O.OlFlo ~ 0.3 pJ. 
Such picojoule latent heat can be readily observed via 
standard ac calorimetry techniques used to measure spe¬ 
cific and latent heats in films: attojoule level results have 
been reported^. We see, therefore, that the latent heat 
Q (1361) associated with the first-order 0 — tt transition in 
SFS trilayer is quite observable. 

We have proposed the general form of the Ginzburg- 
Landau functional to describe SFS trilayer for arbitrary 
phase difference ip between the order parameters in the 
superconducting layers. Calculation of the current-phase 
relation I{(p) shows that the ground state of the SFS 
junction is 0 or tt, and the transition between the 0 and tt 
states appears discontinuous. The current-phase relation 
strongly deviates from the simple sinusoidal one due to 
strong dependence of superconductivity in thin S layers 
on the structure of the pair wave function in the ferro¬ 
magnetic, even at temperatures T near the critical value 
Tc- The characteristic anharmonicity of the current- 
phase relation results in disappearance of the typical non¬ 
monotonic temperature dependence of the critical cur¬ 
rent in a vicinity of the 0 — tt transition. Certainly, the 
anharmonicity of the current-phase relation becomes less 
pronounced for thick ferromagnetic layer, if ^ 

We show that coexisting stable and metastable 0 and 
TT states appear in the vicinity of 0 — tt transition. As a 
consequence, integer and half-integer fluxoid configura¬ 
tions exist in the superconducting loop interrupted by the 
junction. The coexistence of 0 and tt states is manifested 
as two jumps in the dependence of enclosed magnetic flux 
in the loop per period^. 

acknowledgmentsThe authors thank E. Silva, V.V. 
Ryazanov and A.S. Mel’nikov for stimulating discussions. 
This work was supported, in part, by the NanoSC-COST 
(Belgium), Action MP1201 , by ANR (France), Grant 
MASH and by the Russian Foundation for Basic Re¬ 
search. One of the authors (A.V.S.) is supported by 
the Russian Scientific Foundation Grant No 15-12-10020. 

Appendix A: The current—phase relation for rigid 
boundary conditions 

Let us briefly remind of approach developed in Ref. [ij 
to study the problem of the second-harmonic contribu¬ 


tion to the the current-phase relation (l38l) in the frame¬ 
work of the rigid boundary conditions. The general ex¬ 
pression for the supercurrent through a SFS junction is 
given by 


Is = -/o(T/T,) ^ Im {E+F;} , (Al) 

cj>0 


where F+(s, h) = F*{s^ —h), Iq = AnTccSN(0)Df /, S 
is the area of cross section of the junction, and N{0) is 
the electron density of states per one-spin projection. 

The solution of Eqns. m, © for T < Tc and an arbi¬ 
trary phase difference (p is^i 


F{s) ~ a cosh{qs) + b sinh(ps) 


8/c2 


a + — ] [acosh(3fcs) -I- 6sinh(3fcs 


(A2) 

)], 


where the complex wave vectors q and p are different 
from the wave vector k due to nonlinear effects. Just 
below Tc the nonlinear corrections seem to be small and 
wave vectors q and p are determined by the relations 

9 - (a + V4) + ^ (« + 5^4) , (A3) 

P-k-^ia + 5i/4)-^ia + i/A) . (A4) 

The function F+(s) is obtained by replacing b in ex¬ 
pression (IA2I) by —b. Neglecting the influence of the 
E layer on the S layers, one can find the amplitudes a 
and b from the rigid boundary conditions assuming that 
the anomalous Green’s function for the Matsubara fre¬ 
quency uj at the boundary of left (right) S layer coincides 
with the bulk one: Fs{uj,^Sf) = Ab /H, where 

As is the temperature dependent BCS order parameter 
and H -|- lA^p. Using the boundary conditions 


^(±s/) = As(w, ±sf) (A5) 


we get from (IA2I1 , (IA4I1 , (IA4I) the following expressions for 
the amplitudes a and b: 


a = ao + ai, b = bo + bi, 


(A6) 


_ Ab cos(i/?/2) 

H cosh((5/2) ’ 


. Ab sin((/?/2) 
fl sinh(i5/2) ’ 


(A7) 
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Oi = 

bi = 


-ag 


8k^ cosh(5/2) 
-bo 

8k? sinh(5/2) 


[(6q — Aq) (a + 3z/4) (cosh(3<5/2) + 25sinh((5/2)) + i5 (og + 6 q) sinh(i5/2)] 
[(6 q — ag) (a + 3i/4) (sinh(3(5/2) — 25 cosh(^/2)) + id (ag + 6 q) cosh(5/2)] 


(AS) 

(A9) 


where 5 = 5x+i52 and 5x^2 = rf//‘f/i,/ 2 - Substitution of the solutions (IA2I) . (IA6I) . (IA7I) . (IA8I) to the general expression 
for the total Josephson current (lAll) results in the following expression for the amplitude of the first harmonic J1 of 
the current-phase relation 


^ f ^ 

sir. 


Im 


ik 


12 I T, 


j ^i(a + 5j/4) i(a + j/4) (1 — (5/tanh5) I 


sinhJ J 12 V r, y 1 /csinh^ fcsinh^<5 

The second harmonic amplitude is much smaller and described by the following expression 

Iq ^ Ab\^ ^ ( 1 


I 2 = 


(t)M: 


(5 ,a + i/4 ( , ; 5 

2-’ 


192 \ Tc J [ k sinh^ 5 

Actually the Josephson energy (l3^ is the phase dependent contribution to the Ginzburg-Landau energy (|40l) . 


1 / A 


$ 0^2 

P 2 cos 2ip = —-— cos 2ip. 

Attc 


(AlO) 


(All) 


(A12) 


The last equality determines the coefficient (32 in the functional (1401) : 


/32 = - 


192 I T, 


Tco d \ fc sinh^ J [ 2 




sinh J 


sinh J 


(A13) 
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